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Abstract 

Spin orbit interactions are rapidly emerging as the key for enabling efficient current- 
controlled spintronic devices. Much work has focused on the role of spin-orbit cou¬ 
pling at heavy metal/ferromagnet interfaces in generating current-induced spin-orbit 
torques. However, the strong influence of the spin-orbit-derived Dzyaloshinskii- 
Moriya interaction (DMI) on spin textures in these materials is now becoming ap¬ 
parent. Recent reports suggest DMI-stabilized homochiral domain walls (DWs) can 
be driven with high efficiency by spin torque from the spin Hall effect. However, the 
influence of the DMI on the current-induced magnetization switching has not been 
explored nor is yet well-understood, due in part to the difficulty of disentangling spin 
torques and spin textures in nano-sized confined samples. Here we study the mag¬ 
netization reversal of perpendicular magnetized ultrathin dots, and show that the 
switching mechanism is strongly influenced by the DMI, which promotes a univer¬ 
sal chiral non-uniform reversal, even for small samples at the nanoscale. We show 
that ultrafast current-induced and field-induced magnetization switching consists on 
local magnetization reversal with domain wall nucleation followed by its propagation 
along the sample. These findings, not seen in conventional materials, provide essential 
insights for understanding and exploiting chiral magnetism for emerging spintronics 
applications. 

PACS numbers: 75.75.+a, 75.40.Mg, 75.40.Gb, 75.10.Hk 

Keywords: Magnetization Switching, Spin-Orbit Coupling, Spin Hall Effect, Dzyaloshinskii-Moriya interac¬ 
tion 
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Understanding and controlling the current-induced magnetization dynamics in high per¬ 
pendicular magnetocristaline anisotropy heterostructures consisting of a heavy-metal (HM), 
a ferromagnet (FM) and an oxide (HM/FM/O) or asymmetric HM1/FM/HM2 stacks, is 
nowadays the focus of active research 1 -—. Apart from their interest for promising spintronics 
applications, these systems are also attracting growing attention from a fundamental point 
of view due to the rich physics involved in the current-induced magnetization switching 
(CIMS)U- and in the current-induced domain wall motion (CIDWM)- - — . Indeed, the com¬ 
bination of a HM and a thin FM him gives rise to new phenomena which normally vanish 
in bulk, but play an important role as the thickness of the FM is reduced to atomistic size. 

Current-induced torques arising from spin-orbit phenomena can efficiently manipulate 
magnetization. In particular, the Slonczewski-like spin-orbit torque (SL-S()T)-U— can switch 
the magnetization from up (d) to down (d) states and vice versa under the presence of small 
in-plane fields. The SL-SOT is expressed as 

TSL = -70 H SL rh x (m x a) (1) 

where 70 is the gyromagnetic ratio, m the unit vector along the magnetization, cr = u z xuj = 
u y the unit vector along the polarized current which is perpendicular to both the easy axis 
(z) and current direction given by Uj, and Hsl parameterizes the torque. CIMS in ultrathin 
Pt/Co/AlO, where the Co layer is only 0 . 6 nm thick (around three atomic layers), was 
experimentally observed first by Miron and coworkersi, where the switching was attributed 
to SL-SOT due to the Rashba heldr-i^. The Rashba effect would generate both held-like (FL- 
SOT ) 17 ' 18 and Slonczewski-like (SL-SOT)^^ spin-orbit torques. Similar to the conventional 
spin transfer torque (STT)^, both Rashba FL and SL SOTs have magnitudes proportional 
to the spin polarization of the current (P) flowing through the FM, and therefore, they 
are expected to be negligible for an ultrathin FM, as reported in experimental studies^'—. 
Indeed, Liu et al.- studied CIMS in Pt/Co/AlO, similar to the study by Miron et al.- but 
they did not find any significant dominant Rashba FL torque, and therefore the Rashba 
contribution to the SL-SOT should be even vanishingly small. This was also the conclusion 
from switching experiments in asymmetric Pt/Co/Pt^ and for Pt/CoFe/MgO^. Instead of 
the Rashba SL-SOT, the switching is consistent with an alternative SL-SOT based on the 
spin Hall effect (SHE)^^. The SL-SOT due to the SHE is physically distinct from other 
torques STTs and Rashba-SOTs: it is independent of P because it arises from the spin 
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current generated in the HM, rather than the spin polarization of the charge current in the 
FM. 

The key to the existence of the SOTs is a high spin-orbit coupling combined with struc¬ 
tural inversion asymmetry (SIA) in these heterostructures: if the top and bottom inter¬ 
faces/layers sandwiching the FM were completely symmetric, all the mentioned effects should 
cancel out. However, not only the SIA plays a role in these current-induced magnetization 
dynamics but, it can also influence the static magnetization state through the interfacial 
Dzyaloshinskii-Moriya interaction (DMI) 2 - - —. In systems with SIA, the interfacial DMI is 
an anisotropic exchange contribution which directly competes with the exchange interaction, 
and when strong enough, it promotes non-uniform magnetization textures of a definite chiral¬ 
ity such as spin helixes^, chiral domain walls (DWs)^-il and skyrmions^ 2-34 . In particular, 
the experiments on current-induced DW motion along Pt/Co/A10 7 - or Pt/CoFe/MgO^> 41 . can 
be explained by the combined action of the DMI and the SHE. The strong DMI in these Pt 
systems is the responsible of the formation of the Neel walls with a given chirality, which 
are driven by the SHE^ - — . However, the influence of the DMI on the CIMS has not been 
explored nor is yet well-understood, due in part to the difficulty of disentangling spin torques 
and spin textures in nano-size confined dots. 

On the other hand, experiments on CIMS in these asymmetric multilayers are usually 
interpreted in the framework of the single-domain model (SDM) which neglects both the 
exchange and DMI contributions, and only a few recent studies in extended samples at 
the microscale (15/im x 1.2/irn) have considered the non-uniform magnetization by full 3D 
micromagnetic simulations^ - —. Here we focus on CIMS of a ultrathin Pt/Co/AlO with 
in-plane dimensions two orders of magnitude below lOOnm). Although these dimensions 
should be amenable for the uniform magnetization description, our study indicates that the 
DMI is also essential to describe the CIMS at these dimensions, which occurs through chiral 
asymmetric DW nucleation and propagation. We analyze the key ingredients of the switching 
and confirm that a full micromagnetic analysis is necessary to describe and quantify the spin 
Hall angle under realistic conditions. 

RESULTS 

The considered heterostructure here consists on a thin ferromagnetic Co nanosquare with 
a side of L = 90nm and a thickness of L z = 0.6nm sandwiched between a AlO layer and 
on top of a Pt cross Hall (Fig. Oja)). The thickness of the Pt layer is 3nm. Typical high 
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PMA material parameters were adopted in agreement with experimental values^ 1 ^. Details 
about the physical parameters can be found in Methods. 

Cuasi-uniform current-induced magnetization switching in the absence of 
the DMI: single domain approach and micromagnetic results. The current in¬ 
duced magnetization dynamics under static in-plane longitudinal held B = Bu x and cur¬ 
rent pulses j(t) = j{t)u x is studied from both Single Domain Model (SDM) and full mi¬ 
cromagnetic Model (/iM) points of view (see Methods). We first review the CIMS in 
the framework of the SDM, where the magnetization is assumed to be spatially uniform 
( m(t ) = (m x (t ), rn y {t ), m z {t ))). Within this approach the conventional symmetric exchange 
and interfacial DMI are not taken into account (O — 0). In the absence of in-plane fields 
(B = 0) or thermal fluctuations, with the magnetization initially pointing along the easy z- 
axis (m(0) = ±u z , t / |), a moderate current j along the longitudinal direction (x-axis) only 
generates an effective SffE held along the :r-axis which does not promote the out-of-plane 
magnetization reversal ( Hsh = HshU )ffi x a = —HsH(j)m z u x ). ffowever, in the presence 
of a longitudinal held B = Bu x below the saturating in-plane held ( B sat ~ IT), m acquires 
a finite longitudinal component rn x ^ 0 parallel to B , and the current pulse j(t ) generates 
an out-of-plane component effective SHE held H S h,z — Hsh{J])'^x- If jit) is parallel to B 
(either j(t) > 0 and B > 0 as in Fig. [U(c), or j(t ) < 0 and B < 0 as in Fig. [0(e)), and 
their magnitudes are sufficiently strong, the magnetization is stabilized pointing parallel to 
the out-of-plane component of Hsh,z = HsH r m x < 0, i.e. along the — z-axis (Fig. [Qc) and 
(e)). On the contrary, if the held and the current pulse are anti parallel to each other (either 
j(t ) > 0 and B < 0, or j(t) > 0 and B > 0), m is stabilized along the + 2 -axis (Fig. [0(d)). 

Fig. HKf) shows the 3D magnetization trajectories for CIMS starting from the up state 
(t) with j(t) > 0 and B > 0 for 9sh = 0.2 in the absence of DMI (,D = 0). In this case, 
the reversal occurs via quasi-uniform magnetization precession, and therefore, the SDM 
reproduces accurately the magnetization dynamics (solid red line in Fig. [U(f)) computed 
from a full /iM point of view (black dots in Fig. [0(f)), conhrming the validity of the uniform 
magnetization approach in the absence of DMI (0 = 0). 

The SDM stability phase diagrams showing the terminal out-of-plane magnetization di¬ 
rection as function of B and j{t ) (with t R = t F = 200ps, r = 20ps and different amplitudes 
j a ) are depicted in Fig. [ljg) and (h) for a high Osh = 0.2 and a more realistic Osh = 0.11 
value of spin Hall angle respectively. These results were computed at room temperature 
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by averaging over 10 stochastic realizations. The same results were also obtained at zero 
temperature (see open circles in Fig. [Tfg) and (h)). Note that Osh = 0.2 is around twice 
the value experimentally deduced for the Pt/Co from efficiency measurements^, where Osh 
was estimated 0.11. Therefore, these experiments^ cannot be reproduced by the SDM unless 
unrealistic values of Osh ~ 0.4 are assumed^. As it will be shown later, the key ingredient to 
achieve quantitative agreement is the presence of DMI, which can only be taken into account 
in a full fiM analysis. 

Non-uniform magnetization patterns and current induced magnetization 
switching (CIMS) in the presence of finite DMI: micromagnetic results. Al¬ 
though the SDM could qualitatively describe the stability phase diagrams, it fails to provide 
a quantitative description of the experiments^, and the spatial magnetization dependence 
(m(r, t)) needs to be taken into account for a realistic analysis. Indeed, it has been argued 
that the Dzyaloshinskii-Moriya interaction (DMI) arises at the interface between the HM 
(Pt) and the FM (Co) layers^’ 38 . In particular, it was confirmed that apart from SL-SOT 
due to the SHE, also the DMI is a key ingredient in governing the statics and dynamics 
of DWs along ultrathin FM strips sandwiched in asymmetric stacks^’ 30,30 . Similarly to the 
conventional symmetric exchange interaction ( H exc h ) responsible of the ferromagnetic order, 
the interfacial DMI effective field Hdmi is only different from zero if the magnetization is a 
non-uniform continuous vectorial function m = rh(r,t). Apart from promoting non-uniform 
magnetization textures of a definite chirality in the bulk of the FM, the interfacial DMI 
also imposes specific boundary conditions (DMI-BCs) at the surfaces/edges of the sample^. 
Indeed, for finite DMI (D ^ 0), the DMI-BCs ensure that the local magnetization at the 
edges rotates in a plane containing the edge surface normal n, and therefore, in a finite- 
ferromagnetic dot the uniform state is never a solution, so the SDM does no longer apply. 
Further details of the /iM are given in Methods. 

Non-uniform equilibrium states under B = 0 and B ^ 0 in the absence of 
current. In the equilibrium state at rest (B — j — 0), the average magnetization (< m{r) >, 
where < ... > represents the average in the FM volume) points mainly along the easy 
axis, either along +u z (d, Fig. [21(a)) or —u z (d, Fig. [2](b)). However, rh(f) deviates from 
this easy axis direction at the edges (see Fig. I2D- For the up d state (< m >~ +«,), 
the local magnetization m(x,y ) depicts a finite longitudinal component (m x ^ 0), with 
rh(0,y) = +\m x \u x + m z u z and rh(L,y) = —\m x \u x + m z u z at the left (0, z/) and at the right 
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(L,y) laterals respectively (see Fig. [21(a)). Similarly, m(x,y) has a non-zero transversal 
component ( m y ^ 0), with m(x, 0) = + \m y \u y + m z u z and with m(x, L ) = — \m y \u y + rn z u z 
at the bottom (x, 0) and top (x, L ) edges respectively. Instead of pointing inwards (Fig. 
[21(a)), the directions of the in-plane components (■ m x ,m y ) at the edges reverse to outwards 
for the j. state (< m >~ — u z , Fig. [21(b)). The deviations from the perfect out-of-plane state 
are maximum at the edges and decrease over a distance given by ~ ^ toward to the sample 
center. 

A moderate positive longitudinal field B = Bu x well below the in-plane saturating held 
slightly modifies the out-of-plane magnetization in the central part of the FM sample, but 
it introduces significant changes in the local magnetization at the edges, as it can be seen 
in Fig. [21(c)-(d). A finite longitudinal component m x parallel to B arises at both bottom 
and top transverse edges (y : 0, L) (see m y (x,y) in Fig. [21(c)- (d)). Importantly, the effect of 
the positive held B = Bu x with B > 0 is opposite at the longitudinal left (0,y) and right 
(L, y ) edges. Whereas B > 0 supports the positive longitudinal magnetization component 
at the left edge (0 ,y), it acts against the negative longitudinal magnetization component at 
the right edge (L, y) for the t" state, as it is clearly seen in Fig. [21(c). For the state, B > 0 
supports the positive m x (L,y) and acts against the negative m x (0,y) (Fig. [21(d)). 

Non-uniform CIMS from t" to j. with B > 0 and j(t) > 0 for finite DMI (D ^ 0). 
Since for finite DMI (D ^ 0) the equilibrium states of Fig. [21 depict non-uniform magne¬ 
tization patterns m(x^y), and the SHE effective field depends on the local magnetization 
(. Hsh(x,v ) = HsHi'n(x, y) x u y ), the magnetization dynamics must be also non-uniform, 
even for the small nano-sized confined dots with L = 90nm with strong DMI. The non- 
uniform magnetization dynamics under static longitudinal field (B = ±300mT) was studied 
under injection of current pulses j = j(t)u x (Fig. [U(b)) with tn = tp = 200ps, r = 20ps 
and j a = ±3.7 x 10 12 A/m 2 (corresponding to an uniform current / = 1.2mA through the 
Pt/Co section, 3.6 x 90nm 2 ) by yM solving the dynamics equation (Methods). The value 
for the spin Hall angle is 9sh — 0.11 as deduced experimentally by Garello et ah- for similar 
samples. The temporal evolution of the Cartesian magnetization components averaged over 
the volume of the FM (< rrq > (t) with i : x,y,z) and the current pulse temporal profile 
are shown in Fig. [3] for different combinations of B and j a which promote the CIMS 
from j" to (B > 0 and j > 0), and from j, to t (B > 0 and j < 0). Representative transient 
magnetization snapshots during the CIMS are also shown in Fig. [3l which clearly indicate 
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that the switching is non-uniform as opposed to SDM predictions. 

We focus our attention on the CIMS from d to d with B > 0 and j(t ) > 0 (left graphs 
in Fig. [3]) in the presence of strong DMI (D = 1.4mJ/m 2 ). The temporal evolution of the 
Cartesian magnetization components over the ferromagnet volume (< mj > with i : x,y,z) 
is shown in Fig. El a), whereas representative transient magnetization snapshots are shown 
in Fig. El)b)-(f). The reversal takes place in two stages. The first one consists on the 
magnetization reversal at the top left corner of the square resulting in DW nucleation, and 
the second one occurs via current-driven domain wall (DW) propagation from the left to 
right due to the SHE. Apart from the snapshots of Fig. El(b)-(f), these two stages are also 
evident in the temporal evolution of the out-of-plane magnetization < m z > shown in Fig. 
EKa). From t — 0 to t fa 200ps, < rn z > decreases gradually, whereas it decreases almost 
linearly from t zs 200ps to t ~ 300ps, consistent with the current-driven DW propagation 
where its internal structure is seen in Fig. El(d). 

The magnetization reversal during the first stage is non-uniform due to the DMI imposed 
boundary conditions (DMI-BCs, see Methods), but to understand in depth the underlaying 
reasons, it is needed to take into account the chiral-induced non-uniform magnetization 
(m(x, y ) = m x (x, y)u x + m y (x, y)u y + m z (x , y)u z ) in the presence of the applied field (B > 0) 
and current (j >0). As it can be seen in Fig. [2](c) or in Fig. E(b), B > 0 and DMI- 
BCs support the positive longitudinal magnetization component (m x (x,y)) at the left-edge 
(x,y) = (0 ,y), whereas the negative m x (x,y ) is very small at the right edge (x, y) = (L,y). 
An schematic view of the local equilibrium magnetization at relevant locations is shown in 
Fig. EKg) for the d state under B > 0 and zero current. The effective SHE field is also non- 
uniform: H SH (x,y) = H SH ,x(x,y)u x + H SH , z {.x,y)u z with H SH:X (x,y ) = -Hs H (j)m z (x,y) 
and HsH,z(x,y ) = HsH(.j) m x( x ,y)- As the out-of-plane component HsH,z(x,y) is negative 
(note that H S hU ) < 0 for j[t) > 0) and proportional to the local 7n x (x,y), which is 
maximum and positive at the left edge (( x,y ) = (0,r/)), the reversal starts from the left 
edge (see Fig. El(h)). However, in addition to this asymmetry along the longitudinal x-axis 
imposed by the DMI-BCs and supported by B (left vs right edges), other chiral asymmetry 
arises along the transverse y- axis in the left edge: the reversal is first triggered from the 
top left corner ((x,y) = (0,L)), whereas the local CIMS is delayed at the bottom-left 
corner ((x, y) = (0,L)), as it clearly seen in Fig. El(c). The reason for this transverse 
asymmetry relies in the different direction of local torque at the initial state (Fig. Eli))- 


8 


The relevant torque is the one experienced by the local magnetization at the left edge 
rh(0,y) due to Hsh(x,i/), which is also supported by B = y 0 H x u x : r(x,y ) = —70 rh(x,y) x 
(Hsh(x,v) + H x ) = r x u x + Tyily + r z u z . As the local transverse magnetization jn y (x,y ) 
has different sign at the top (m y (0,L) < 0) and bottom (m. y (0,0) > 0) corners of the left 
edge, both the longitudinal component (r x (x, y) = —70 m y (x, y)m x (x, y)H S H ) and the out-of- 
plane component of this torque ( r~(x,y ) = —70 (m y (x,y)m z (x,y)HsH — m y( x ,y)H x )) point 
in opposite directions at the top and the bottom corners of the left edge (see Fig. my 
The relevant component of t(x, y) to understand the local reversal is the out-of-plane one: 
as r 2 (0, L) < 0 at the top left corner but r 2 (0, 0) > 0 at the bottom left corner, the reversal 
is firstly triggered from the top corner, where t z (0,L) opposes to the initial out-of-plane 
component of the magnetization (d). Once the local reversal is achieved at the top left 
corner, the switching expands from left to right and from top to bottom: the local in-plane 
magnetization at the bottom left edge rotates clockwise due to Hg H , z (0 ,0) < 0, and once 
777^(0,0) becomes negative, also r 2 (0,0) < 0 promotes the local reversal. 

When all points at the left edge have reversed their initial out-of-plane magnetization 
(m z (0,y) < 0) a left-handed (D > 0) down-up DW emerges, separating the reversed (with 
d) from the non-reversed (with d) zones. Note that once the local magnetization has reversed 
its initial out-of-plane direction, it experiences little torque due to Hsh (see Supplementary 
Information), so it is stable for the rest of the switching process, which takes place by 
current-driven DW propagation during the second stage. 

The internal structure of the propagating DW is shown in Fig. [3](d). Even in the presence 
of the longitudinal field (B > 0), its internal moment [mow) and its normal (now) do not 
point along the positive x-axis, and the DW depicts tilting or a rotation of its normal due 
to the SHE current-driven propagation. The DW tilting has been experimentally observed 
in the absence of in-plane field under high currents^, and theoretically studied, both in the 
absence and in the presence of in-plane fields, in elongated strips along the axaxis^^h—. If 
the only driving force on the down-up DW (dd) were a strong positive (negative) current 
j > 0 (j < 0) with B — 0, both mow and now would rotate clockwise (counter-clock 
wise)^. Here, we observe that the DW tilting is also assisted during the DW nucleation due 
to the DMI-BCs, B and Hsh- B > 0 would support the internal longitudinal magnetization 
of the left-handed down-up DW if its normal points along the x-axis (now — +u x ), as it 
would be the case of current-driven DW motion along an elongated strip along the x-axisA^. 
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However, due to the non-uniform local CIMS at the left edge in our confined dots, the DW 
normal has a non-zero negative transverse component (riow.y < 0) for B > 0 and j > 0. 
As it is shown in the 270ps-snapshot of Fig. [31(d) , in addition to a positive longitudinal 
component (mjjw,x > 0), the internal DW moment also has a no-null negative transverse 
component ( mr>w,y < 0). Note that the direction of both n DW and rh DW during the DW 
propagation is also the direction of the local magnetization at the top-left corner, where the 
reversal was initially launched (see Fig. Hc),(i)). 

The full magnetization switching is completed before the current pulse has been switched 
off (see Fig. Ela)), when the propagating down-up DW (^t) reaches the right edge. Due to 
the DW tilting, the reversal occurs first at the top right corner ((ay y ) = (L, L )) with respect 
to the bottom right corner ((x,y) = (L, 0)) (see 320ps-snapshot of Fig. E](e)). Although this 
second stage, consisting on current-driven DW propagation, is similar to the one already 
explained for elongated thin strips as driven by the SHEAhiZd^ the DW nucleation during 
the first stage has not been addressed so far for such small nano-sized confined dots, and 
as it was explained above it is mainly due to the longitudinal held B which supports the 
longitudinal magnetization component at the left edge imposed by the DMI-BCs. 
DISCUSSION 

Universal chiral promoted current-induced magnetization switching (CIMS) 
in strong DMI systems. The CIMS from to J, can also be achieved if both B and 
j reverse their directions (B < 0 and j < 0). As it is straightforwardly understood from 
the former description, in this case the reversal is triggered from the bottom right corner 
(( x , y) = (L, 0), where t 2 (L, 0) < 0 opposes to the initial f out-of-plane magnetization), and 
an up-down DW is driven toward the left (not shown). The CIMS from J, to t under anti 
parallel held B > 0 and current j < 0 is shown at the right panel of Fig. [3](j)-(r). 

In general, the CIMS can be described as follows: (i) the initial out-of-plane magnetization 
direction (f or ),) determines the direction (inwards or outwards) of the local in-plane m at 
the edges imposed by the DMI-BCs. (ii) The longitudinal held B supports the longitudinal 
in-plane magnetization component ( m x ) at one of the two lateral edges, and acts against it 
at the opposite one. (iii) For the favored lateral edge, the local magnetization reversal is 
triggered at the corner where the out-of-plane torque t z due to Hsh and B opposes to the 
initial out-of-plane magnetization component (m z ). After that, the reversal also takes place 
in the middle part of the selected edge, and finally, the other corner is also dragged into 
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the reversed region with the formation of a tilted DW. (iv) The CIMS is completed by the 
current-driven DW propagation. 

Also remarkable is the fact that for the same current pulses as in Fig. [3] the CIMS is 
not achieved in the framework of the SDM if a realistic value for the spin Hall angle is 
adopted ( 9sh = +0.11}£, and the same limitation was also observed by full fiM simulations 
in the absence of the DMI (D = 0). All these simulations point out that, even for the 
small confined dots considered here (L = 90nm), the strong DMI and the BCs imposed 
by it are essential to describe the CIMS driven by the SHE from both quantitative and 
qualitative points of view. The DMI-triggered switching (D ^ 0) was also studied for 
other ultrathin ( L z = O.611111) squares ( L x = L y = L) with different in-plane dimensions 
(lOnm < L < 300nm) and reversal mechanism remains similar to the one already described 
and depicted in Figj3l Note that the smallest evaluated side (L = lOnm) is small than 
the minimum side required to achieve thermal stability ( L min m 25nm) according to the 
conventional criterion: in order to maintain sufficient stability of the data storage over at 
least five years, the effective energy barrier given by Eb = KV (with K fa 4.35 x 10 5 J/m 3 the 
effective uniaxial anisotropy constant from Ref.-, and V = E 2 L Z the volume of the sample) 
should be larger than 55K B T, where and K B Boltzmann constant. The reversal was also 
similar under realistic conditions including disorder due to the edge roughness and thermal 
effects (see Supplementary Information). Moreover, this chiral CIMS, either from "f to or 
from j. to d, does not change when the FM Co layer is patterned with a disk shape (see 
Supplementary Information). It was also verified that this non-uniform reversal mechanism, 
consisting on DW nucleation and propagation, does not depend on the specific temporal 
profile of the applied pulse, provided its magnitude ( j a ) and duration (r) are sufficient to 
promote the complete reversal for each L. 

Chiral nature of the field-induced magnetization switching (FIMS). An anal¬ 
ogous CIMS mechanism to the one described here for nano-size samples (~ 90nm x 
90nm) was recently observed by Yu et al. 49 using Kerr microscopy for an extended 
Ta(5nm)/CoFeB(lnm)/TaO(1.5nm) stack with micro-size in-plane dimensions (20/im x 
130/mi). In that work, right-handed DWs (D < 0) were nucleated assisted by the in-plane 
field and displaced along the current direction due to the negative spin Hall angle of the Ta. 
More recently, Pizzini et ah— also used Kerr microscopy to visualize the asymmetric chiral 
DW nucleation under in-plane field and its subsequent propagation along extended (~ 70/mi) 
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Pt(3nm)/Co(0.6nm)/A10(2nm) thin-films driven by out-of-plane field ( B oop = B z u z ). Simi¬ 
lar to our study, starting from the up state (t), a positive (negative) in-plane field (B = Bu x ) 
promotes the local magnetization reversal at the left (right) edge, which was propagated to 
the right (left) driven by a negative out-of-plane held B z < 0. Their images indicate the 
nucleated DW has a left-handed chirality and it propagates without significant tilting due 
to the extended unconfined in-plane dimensions (L y ps 70/nn). In order to understand these 
observations, the field-induced magnetization switching (FIMS) has been also studied for 
confined small squares with L = 90nm (the same geometry as in the former CIMS analysis) 
and others with lateral dimensions one order of magnitude larger (L = lOOOnm). Static 
longitudinal fields B = Bu x with B > 0 and B < 0 (\B\ = 300mT) are applied along with 
short out-of-plane held pulses with B oop (t ) = B z (t)u z with \B Z \ = 310mT, t r = tp = 200ps 
and r = 20ps (the temporal prohle of this pulse is the same as for the current-induced mag¬ 
netization switching). The results for the confined L = 90nm square dot are shown in Fig. 
[4] for different combinations of the initial state (t and J,), in-plane static held B (B > 0 and 
B < 0) and out-of-plane held pulse ( B z {t ) > 0 and B z (t ) < 0). Similarly to the CIMS, the 
FIMS starts from an edge selected by the direction of H, with an even more evident chiral 
asymmetry between the two corners. Note again that the corner where the reversal starts 
has a transverse magnetization component (m y (x,y)) pointing in the same direction as the 
transverse internal magnetization of the nucleated DW. Once the local switching has been 
triggered, the reverse domain (pointing along the opposite 2 direction with respect to the 
initial state) expands asymmetrically along the longitudinal (a:) and transverse ( y ) directions 
(see for instance snapshots at t = 270ps and t = 320ps in Fig. [4]). Although here just a 
quarter-of-bubble is developed due to the confined shape at the corner, this asymmetric field- 
driven chiral expansion is similar to the one recently observedr^*^ in extended thin films. 
Moreover, our study also points out a qualitative difference between the current-driven and 
the field-driven nucleation: while the first one is driven by a non-uniform SHE out-of-plane 
effective held ( B S H,z(x,y ) = /i 0 H S H,z(x,y) which depends on local m x (x,y)), the second one 
is promoted by a uniform out-of-plane held B z . Therefore, the current-induced nucleated 
DW propagates along the current direction (a:-axis, see yellow arrows in Fig. |3](d) and (m)), 
whereas the field-driven DW expands radially from the corner (see yellow arrows in Fig. 
HD- Nevertheless, the fact that similar chiral local magnetization reversal occurs also at the 
corners of nano-size confined dots (~ lOOnm and below) clearly confirms the universality of 
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the chiral reversal mechanism in these nano-size confined dots with strong DMI. 

On the other hand, the Kerr images by Pizzini et ah— do not show the corners of their ex¬ 
tended thin-hlm (which is unconfined along the transverse y-axis) which are precisely where 
our modeling points out additional chiral asymmetry in the DW nucleation for confined 
dots (Fig. [3J. Moreover, in their thin-films the field-driven DW does not depict tilting. In 
order to contrast these observations with our yM predictions, the field-driven nucleation 
and propagation in an confined square dot has been also analyzed here, but with lateral 
in-plane dimensions one order of magnitude larger (L = lOOOnm). We note that as L is 
increased to the microscale, the nucleated DW is almost straight, with its normal oriented 
along the x-axis (no DW tilting), in the middle part of the nucleating edge (far form the 
corners). However, an asymmetry between the top and bottom corners is still present even 
for L = lOOOnm (see Supplementary Information): the reversal from "f to J, (from j. to d) is 
anticipated at the top-left (top-right) corner with respect to the bottom one under B > 0 
and B z < 0 (B < 0 and B z < 0). This chiral asymmetry at the corners of the extended 
micro-size sample is similar to the observed for a confined dot (see. Fig. SD, and although 
it has not been addressed before, it could be observed by high resolution techniques^. 

CIMS in confined nanodots with rectangular shape. The CIMS was also studied in 
rectangles with different in-plane aspect-ratios 1 < L x /L y < 4 (Fig. [5](a)-(c)). The thickness 
is fixed ( L z = 0.6nm) as before. Again the switching takes place by DW nucleation followed 
by its current-driven propagation along the x-axis, which further supports the universality of 
the reversal mechanism in systems with strong DMI. In this case, the nucleation takes place 
during the first 200ps independently of the rectangle aspect-ratio L x /L y , but the critical 
pulse duration (r) for fixed j a and t R = tp, increases linearly with L x /L y (see the inset in 
Fig. [5](c)), an prediction which could be experimentally validated to estimate both the spin 
Hall angle (Osh) and the DMI parameter ( D ) if the rest of material parameters (M s , A, K u , 
a) are known by other means. 

Comparison to experiments of current-induced magnetization switching. Al¬ 
though our study goes further than a mere comparison to available experimental results, it is 
interesting to show how the non-uniform CIMS can explain quantitatively the experimental 
measurements by considering realistic material parameters (see Methods and Supplementary 
Information). With the aim of providing an explanation of experimental observations^ for 
the ultrahin Co square with L = 90nm in a Pt(3nm)/Co(0.6nm)/A10(2nm) stack, we have 
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repeated the former study for several values of the applied field ( B) and different different 
magnitudes of the current pulse (j a ). The rise and fall times (£r = tp — 200ps) and the 
duration (r = 20ps) of the pulse were maintained fixed as in the experimental study^. Here 
we consider the up state (< m z >~ +1, d) as the initial one. For each ( B,j a ), the switching 
probability at room temperature was computed as the averaged over 10 stochastic realiza¬ 
tions. Realistic conditions were taken into account by considering random edge roughness 
with characteristic sizes ranging from 0.5nm-5nm (see Methods). The /jM results are col¬ 
lected in Fig. [6] which indicates a good quantitative agreement with recent experimental 
measurements^. 

It was verified that the CIMS mechanism (local magnetization reversal with DW nucle- 
ation and subsequent current-driven propagation) remains qualitatively unchanged even un¬ 
der these realistic conditions (see Supplementary Information). Moreover, although marginal 
discrepancies between these fiM data (Fig. [6]) and the experimental results shown in Fig. 
2(d) of ref. 6 ' can be seen, the quantitative agreement is remarkable considering similar ma¬ 
terial parameters as inferred experimentally^: M s = 8.7 x 10 5 A/m, A = 1.6 x 10 -11 J/m, 
K u = 8 x 10 5 J/m 3 , a = 0.3, Osh = 0.11 and D = 1.4mJ/m 2 (see Supplementary Information 
for detailed justification of these inputs). Note that with the SDM a quantitative agreement 
with the experimental data was only achieved with unrealistic values of the ( Osh = 0.4)&. 
Note that the DMI parameter D = 1.4mJ/m 2 was not determined experimentally^, but the 
fact that this value D = 1.4mJ/m 2 provides reasonable quantitative agreement with their 
experiments, and that this value is also in good quantitative agreement with very recent es¬ 
timations by other means for similar Pt/Co/AlO systems 7 constitute additional evidences 
that our modeling is compatible with the dominant physics behind these CIMS processes. 

Conclusions. In summary, the current-driven magnetization switching in ultrathin 
HM/FM/Oxide heterostructures with high PMA and strong DMI has been studied by means 
of full micromagnetic simulations. Even for the small in-plane dimensions (~ lOOnm), the 
analysis points out that the magnetization reversal mechanism is non-uniform. It starts by 
local magnetization reversal induced by the SHE and assisted by the in-plane held in col¬ 
laboration with the DMI boundary conditions. The longitudinal held and the DMI imposed 
boundary conditions select the lateral/edge and the specihc corner at which the nucleation 
is triggered, where the relevant torques due to the SHE and the longitudinal held accelerate 
the local reversal. After that, the switching is completed by current-driven domain wall 
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propagation driven by the SHE, where the current direction determines the direction of the 
wall motion, and the internal magnetization of the propagating wall points closely to the 
local magnetization at the selected corner where the reversal was initially launched. Similar 
nucleation and propagation mechanisms were also observed under out-of-plane fields, con¬ 
firming again the chiral-triggered magnetization reversal in these nano-size confined dots. 
These results clearly exclude the single domain approach as a proper model to describe 
these switching experiments, and therefore, the estimations of the spin Hall angle based 
in this oversimplified model should be revised by adopting a much more realistic full 3D 
micromagnetic approach. Moreover, by analyzing the switching under realistic conditions 
including disorder and thermal effects, it was found that the mechanism is universal, and for 
instance, it could be used to the quantify both the DMI and the spin Hall angle by studying 
the reversal of ferromagnetic layers with different length for fixed width and thickness. As 
the reversal mechanism occurs in a reliable and efficient way, and more importantly, as it is 
also highly insensitive to defects and thermal fluctuations, our results are also very relevant 
for technological recording applications combining non-volatility, high stability, ultra-dense 
storage and ultrafast writing. 

METHODS 

Magnetization dynamics under SOT due to the SHE. Under injection of a spatially 
uniform current density pulse along the x-axis j(t) = j(t)u x (see its temporal profile in Fig. 
QJb)), the magnetization dynamics is governed by the augmented Landau-Lifshitz Gilbert 
eq. 


dm 

dt 


~ 7o H SH rh x (m x a) 


( 2 ) 


where m(f,t) = M(r,t)/M s is the normalized local magnetization with M s saturation mag¬ 
netization, 70 is the gyromagnetic ratio and H e ff is effective held derived from the energy 
density of system ( H e ff — — m:)- The term in equation (J2J) represents the preces- 

sional torque of m around H e ff + H t h, where H t h is the thermal held representing the effect 
of thermal huctuations at finite temperature. H t h is a white-noise Gaussian-distributed 
stochastic random process with zero mean value (its statistical properties are given below). 
The second term in equation ([2]) is the damping torque with a the dimensionless Gilbert 
damping parameter. The last term in equation fl2]) is the SL-SOT from the spin Hall effect 
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(SHE), where a = u y is the unit vector pointing along the direction of spin current polar¬ 
ization due to the SHE in the Pt layer, and Hsh represents the magnitude of the effective 
spin Hall held Hsh = H s n'm x a given by 


H SH it) = 


hQsH.j(t) 


(3) 


2 en 0 M s L z 

where L z is thickness of the FM layer, Ti is Planck’s constant, e < 0 is the electron charge 
and j(t) is the instantaneous value of the electrical density current. As in the experiment by 
Garello et al.—, the current is assumed to how uniformly through the HM/FM bilayer (see 
Supplementary Information for additional discussion). Osh is the Spin Hall angle, which is 
dehned as the ratio between the spin and charge current densities. 

Single Domain Model (SDM). If the magnetization is assumed to be spatially uniform 
{fh(t) = m x (t)u x + m y (t)u y + m z {t)u z ), the deterministic effective H e ff held in equation ([2]) 
only includes the PMA anisotropy, magnetostatic and Zeeman contributions H e ff = Hpma + 
Hdmg + H ext . The Zeeman contribution due to the longitudinal held is H ext = ( B//i 0 )u x . 
The uniaxial PMA anisotropy effective held is 


f} 2A u 

Hpma — — 

IM)M s 


and the demagnetizing held in the SDM approach is expressed as 


(4) 


Hdmg HI S NHI HI s ( N x Hl x Tl x NyHlylly T N z Hl z U z ) (5) 

where N is the diagonal magnetostatic tensor with N x = N y = 0.011 and N z = 0.975 being 
the self-magnetostatic factors^ for L = 90nm and L z = 0.6nm. 

The thermal held H t h(t ) is a stochastic vector process whose magnitude is related to the 
temperature T via the fluctuation-dissipation theorem 41 . 


Hthit) 




2<aKbT 
y,/M, M S V At 


( 6 ) 


where Kp is the Boltzmann constant, V = L x L y L z is the volume of the sample, At is the 
time step, and ff(t) = (r] x , rj y , rj z ) is a Gaussian distributed white-noise stochastic vector with 
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zero mean value (< rji(t) >= 0 for i : x,y,z ) and uncorrelated in time (< 'rji(t)rjj(t') >= 
SijS(t — t'), where 8,^ is the Kronecker delta and S(t — f) the Dirac delta). Here < ... > 
means the statistical average over different stochastic realizations of the stochastic process. 
Equation ((2]) was numerically solved with a 4 t,l -order Runge-Kutta scheme with a time step 
of O.lps. 

Micromagnetic Model (/_iM ). When the spatial dependence of the magnetization is 
taken into account ( rh(r,t )), the deterministic effective field H e ff(r,t ) in equation © in¬ 
cludes the space-dependent exchange H exc h(r,t) = ^ (V 2 m) with A the exchange con¬ 

stant, and the interfacial DMI H DMI (r,t) = —^ [(V • rh)u z — Vm 2 ^ where D is 
a parameter describing the DMI magnitude. Both the local Zeeman and PMA uniax¬ 
ial contributions to H e ff(r,t) are computed similarly as in the SDM ( H ext = j^u x and 
Hpma{Ti t) — Note also that in the /j,M the magnetostatic field Hd mg (f,t) 

is also space-dependent on rn(r,t) everywhere. The Oersted field due to the current was 
also taken into account but it was found irrelevant and very small as compared to the other 
dominant contributions in H e ff. (see- 42 ^ for the numerical details). 

In the absence of DMI (D — 0), the symmetric exchange interaction imposes boundary 
conditions (BCs) at the surfaces of the sample 44 so that m{f) does not change along the 
surface ( dfh/dn = 0, where d/dn indicates the derivative in the outside direction normal to 
the surface of the sample). However, in the presence of the interfacial DMI (h ^ 0), these 
BCs have to be replaced by^ 1 ^ 


dm D _ 

_ = -- m X („ X 

where n represents the local unit vector normal to each sample surface. 

In the fiM the thermal held H t h(r,t ) is also a stochastic vector process given by 


(7) 




2aK i >T 


( 8 ) 


joH 0 M s AV At 

where now AH = AxAyAz is the volume of each computational cell and ff(f, t) = (■ rj x , r) y , rj z ) 
is a white-noise Gaussian distributed stochastic vector with zero mean value (< rji(r, t) >= 0 
for i : x,y, z ) and uncorrelated both in time and in space (< t)r}j(r', t') >= S itj 8(t - 

t')5(f—r')). Most of the simulations for perfect samples were performed with a 2D discretiza¬ 
tion using cells of Ax = Ay = 2.5nm in side, and thickness equal to the ferromagnetic layer 
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(L z = O. 611111 ). Several tests were performed with cell sizes of 0.5nm to confirm the numerical 
validity of the presented results. Realistic samples were also studied by considering edge 
roughness using cell sizes of 0.5nm. These realistic conditions are introduced by randomly 
generating edge roughness patterns with different characteristic sizes 0.5nm < D g < 5nm 
at all edges. Equation (J2]) was numerically solved with a 6 ^-order Runge-Kutta scheme 
with a time step of O.Olps by using GPMagnet4^, a commercial parallelized finite-difference 
micromagnetic solver—. 

Material parameters. Typical high PMA material parameters were adopted for the 
results collected in the main text in agreement with experimental values for Pt/Co/AIQ^ 1 ^: 
saturation magnetization M s = 8.7 x 10 5 A/m, exchange constant A = 1.6 x 10 _ 11 J/m, 
uniaxial anisotropy constant K u = 8.7 x 10 5 J/m 3 4 . The spin Hall angle is assumed to be 
Osh = 0.11, also according to experiments by Garello et al.— >&. Note that this value is also 
in the middle of the experimental bounds 0.056 < Osh A 0.16 estimated by Liu et al.— and 
Garello et al.-, and very close to the one deduced irt^. A DMI parameter of D = 1.4mJ/m 2 
is assumed, which is similar to the one experimentally deduced by Emori et ah— . The 
Gilbert damping is a = 0.3 as measured in22. Several tests were also performed by varying 
these inputs within the range available in the experimental literature (see Supplementary 
Information). 
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FIG. 1: Current induced magnetization switching in the absence of DMI (D = 0) | (a) 

Schematic representation of the analyzed heterostructure with the Co layer in blue, (b) Temporal 
variation of the density current pulses j(t ) with rising (fi?), falling (tp) and duration (r) times. 
(c)-(e) Out-of-plane component of spin Hall effective field Hsh,z as a function of the applied 
field B and density current j(t) directions, (f) Magnetization trajectories starting from up state 
(t) to down state ( 4 .) under a static field of B = 300mT and a pulse with tp = tp = 200ps, 
r = 20ps and I = 1.2mA (j a = 3.7 x 10 12 A/m 2 flowing through both the Pt and Co layers) for 
Osh = 0.2. Solid red line depicts Single Domain Model (SDM) results whereas solid black dots 
correspond to full micromagnetic (pM) simulations in the absence of DMI (D=0) for the averaged 
magnetization components over the sample volume (< ... >). (g)-(h) Stability phase diagrams 
indicating the terminal out-of-plane magnetization component m z as a function of B and j(t) for 
tR = tF = 200ps, r = 20ps, and T = 300K as computed with the SDM with a high Osh = 0.2 
(g) and with realistic Osh = 0.11 (h). Open circles denote the transition between switching and 
no-switching at zero temperature. 
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FIG. 2: Non-uniform equilibrium magnetization patterns in the presence of finite DMI 

(D = 1.4mJm -2 ) | Magnetization snapshots depict the deviations of the local magnetization 
rh(x,y) from the perfect out-of-plane direction as due to the DMI-BCs (equation (j7|l) at rest 
(B = j = 0) in the presence of interfacial DMI (D = 1.4mJ/m 2 ) for an state mainly up magnetized 
(< fh >« +u z , t) (a), and for an state mainly down magnetized (< m >« — u z , |) (b). Density 
plots of the longitudinal m x (x,y), transverse m y (x,y) and out-of-plane m z (x,y) configuration are 
shown from top to bottom respectively. Arrows show m z (x,y). (c) and (d) show the equilibrium 
state under a positive longitudinal field B = +300nm for the t and 4- states respectively. 
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FIG. 3: Non-uniform current-induced magnetization switching (CIMS) in the presence 

of DMI (D = 1.4mJm -2 ) | Graphs at the left panel correspond to the to \. switching under 
a positive current pulse (j > 0). (a) Temporal evolution of the Cartesian components of the 

magnetization averaged over volume sample (< m >) for B > 0 and j a > 0. The applied pulse 
j(t) is also shown, (b)-(f) Magnetization m(x,y) snapshots during the f to 4, CIMS. Green box 
in (b) indicates the corner where the switching is triggered as explained in the text and in the 
schemes (g)-(i): (g) shows m(x,y) at different points of relevance to understand the CIMS. Dotted 
arrows indicate the in-plane components of the equilibrium m for B = j = 0, whereas solid vectors 
indicate the equilibrium state under B > 0. B supports the in-plane longitudinal component at 
the left edge (x,y) = (0, y). (h) Scheme of the out-of-plane component ( Hsh,z , hi red) and the 
in-plane longitudinal component {Hsh,x, hi purple) of the SHE effective field ( Hsh(xiV )) at the 
left edge corresponding to (b) and (g). (i) Cartesian components of the local torque (t x . r y . t z ) due 
to B = Bu x and HsH(x,y ) at the relevant left edge: the CIMS is triggered at the top left corner, 
where t z < 0 is opposed to the initial out-of-plane up magnetization. Graphs at the right panel 
(j)-(r) correspond to the 4. to t CIMS under B 25 O but j < 0. Yellow boxes in (d) and (m) indicate 
the internal structure of the current-driven domain wall motion due to the SHE. 
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FIG. 4: Field-induced magnetization switching in a ultrathin (L~ = 0.6nm) square dot 
with L = 90nm | Transient snapshots magnetization m(f) in the presence of interfacial DMI 
(D = 1.4mJ/m 2 ) during the reversal for different combinations of B and B z with \B\ = 300mT, 
\B Z \ = 310mT, tR = tp = 200ps and r = 20ps. The up (|) to down (|) switching is shown for 
(B > 0 ,B Z < 0) and (B < 0 ,B Z < 0) in (a) and (c) panels respectively, whereas the down (4.) to 
up (t) is shown in (b) and (d) for (B > 0 ,B Z > 0) and (B < 0 ,B Z > 0). Green boxes indicate 
the region where the DW nucleation starts for each combination of initial state ('f or 4 ,), B and j. 
Yellow circles indicate the field-driven propagating DW (yellow arrow indicate the direction of the 
reversed domain expansion). Note that the in-plane components of the nucleation region (green 
box) point along close to the internal DW moment (yellow circle) during its propagation. 
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FIG. 5: Current-induced magnetization switching along thin rectangles | CIMS in rect¬ 
angles with L y = 90nm and different aspect-ratio L x /L y : 1, 2, 3,4. The applied field B = +300mT 
and the current pulses have j a = 3.7 x 10 12 A/m 2 and = tp = 200ps fixed, and different dura¬ 
tions r depending on L x /L y (a), (b) Temporal evolution of the out-of-plane component < m z > 
for different rectangles under the pulses shown in (b). Snapshots of the magnetization state at 
t ~ 250ps (c) (DW nucleation) and at t ss 350ps (d). The inset in (b) shows the critical threshold 
for r as a function of the L x /L y . 
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FIG. 6: Quantitative description of experimental results | Micromagnetically computed 
switching probability as function of the applied field B and the current pulse j(t). Similar current 
pulses as in the experiments by Garello et al. 6 are applied: tn = tp = 200ps and r = 20ps are 
fixed, and different magnitudes j a = //(3.6 x 90nm 2 ) are studied (/ = 1.2mA corresponds to 
j a = 3.7 x 10 12 A/m 2 ). Results were computed at room temperature T = 300K by averaging over 
10 stochastic realizations, (a) Switching probability as a function of B for pulses with several 
magnitudes expressed in term of I as in the experimental study. The switching probability is as 
function B and / is depicted by density plot in (b). Dashed white curve in (b) represents the 
threshold between not-switching and switching computed at zero temperature. 
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